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A semiclassical distorted wave (SCDW) model with Wigner transform of one-body density 
matrix is presented for multistep direct {p,p'x) reactions to the continuum. The model uses 
Wigner distribution functions obtained in methods which include nucleon-nucleon correlations 
to a different extent, as well as Woods-Saxon (WS) single-particle wave function. The higher 
momentum components of target nucleons that play a crucial role in reproducing the high- 
energy part of the backward proton spectra are properly taken into account. This SCDW model 
is applied to analyses of multistep direct processes in ^^C(p,p'x), '^^Ca.{p,p' x) and ^'^Zr(p,p'x) in 
the incident energy range of 150-392 MeV. The double differential cross sections are calculated 
up to three-step processes. The calculated angular distributions are in good agreement with 
the experimental data, in particular at backward angles where the previous SCDW calculations 
with the WS single-particle wave function showed large underestimation. It is found that the 
result with the Wigner distribution function based on the coherent density fluctuation model 
provides overall better agreement with the experimental data over the whole emission energies. 



1 Introduction 



Preequilibrium processes in nuclear reactions at intermediate energies are known to be dominated 
by multistep direct (MSD) processes. Among the models proposed in the years, the semiclassical 
distorted wave model [1, 2, 3] has proved its efficiency to describe MSD, especially for reactions 
at intermediate energies. It is based on the DWBA expansion of T-matrix elements and the 
cross section formula has no free adjustable parameter allowing a simple intuitive interpretation. 
The SCDW model has been applied to analyses of inclusive cross sections and spin observables 
in {p,p'x) and {p,nx) reactions over incident energy region of 60 to 400 MeV and wide range 
of mass number [4, 5, 6, 7, 8]. Recently, several extensions and improvements of the SCDW 
model have been made. They include calculations up to three-step processes [4], incorporation 
of a single-particle (s.p.) wave function of target nuclei by using the Wigner transform of a 
one-body density matrix (OBDM) [5], introduction of phenomenological effective mass m* of a 
nucleon in the target nucleus [6] and analyses of spin observables [7, 8] which contain important 
information on an effective interaction in nuclear medium. As a result, SCDW model calculations 
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can well reproduce the experimental data of double differential inclusive cross sections (DDX's) 
for {p,p'x) and {p,nx) reactions. The model, however, underestimated the cross sections at very 
backward angles. This suggests that improvement of the SCDW model might be achieved by 
the inclusion of the high momentum components of the target nucleons. 

The problem of backward proton production in proton-nucleus collisions has been a subject 
of much interest for several years. The main difficulty arises from the fact that the reaction 
involves a very large momentum transfer with only a small energy transfer. One of the models 
which have been applied to solve the problem is based on the correlated two-nucleon-cluster 
picture [9]. A good description of the energy spectra even at the very high energy backward 
protons is achieved with this model. The main reason is that a phenomenological momentum 
distribution with high momentum components was employed [10]. These components were 
found to be quite sensitive to the problem of backward proton production and were considered 
to be an effect due to many-nucleon correlations. Within the model [9] single scattering and 
two-nucleon cluster mechanisms at backward angles have been considered. Although there are 
some similarities between these contributions and corresponding one- and two-step processes, 
the advantage of the SCDW model is that the contribution from three-step process can be 
also obtained. Moreover, in SCDW model the distortion effects for the projectile and observed 
protons are carefully considered within DWBA while the calculated single scattering mechanism 
in [9] has evaluated the distortion effect approximately by the Glauber theory. 

In Ref. [5] the SCDW model was modified so that realistic s.p. wave functions in a finite 
range potential can be used in terms of the Wigner transform of OBDM. Single-particle models 
with harmonic oscillator and Woods-Saxon potentials were used to describe the nuclear states 
for ^^Zx{p,p'x) reactions at 80 and 160 MeV instead of the local density Fermi gas (LFG) 
model used in previous analyses [4, 7]. It was concluded in [5] that the SCDW model with 
realistic s.p. wave functions gives larger one-step cross sections at backward angles, which 
result in better agreement with the experimental data over a wider angular range than the 
model with LFG. The calculations still somewhat underestimated the cross section at very large 
angles. Therefore, apart from the necessity of involving contributions of higher MSD processes 
the inclusion of higher momentum components of target nucleons becomes apparent. In the 
present paper we extend the SCDW model in terms of the Wigner transform of OBDM by 
considering the high-momentum components of nucleon momentum distributions deduced from 
several theoretical methods. This extended SCDW model is applied to ^^C(p,p'x) reaction at 
150 McV [11, 12], 40Ca(p,p'x) reaction at 186 [13] and 392 MeV [14] and 9°Zr(p,p'x) at 160 
MeV [15] and comparison with the experimental data is made. 

The aim of our work is to analyze the angular distributions of the backward scattered protons 
through testing of different models for the nuclear states which account for nucleon- nucleon (NN) 
correlations. Such a systematic study could lead to a further refinement of the SCDW model 
and to a better agreement of the DDX with the experimental data at backward angles due to a 
proper inclusion of high-momentum components of the nuclear Fermi motion. 

The paper is organized as follows. A brief formulation of the SCDW model and the proposed 
WT within different theoretical models are given in Sections II and III, respectively. The results 
of the calculations are presented and discussed in Section IV. The summary of the present work 
is given in Section V. 

2 Outline of the SCDW model 

The formulation of the SCDW model used in the present analysis has been described in details 
elsewhere [5, 7, 8], therefore only the outline is mentioned here. 

Let us consider the one-step process in {N^ N'x) reactions in which a target nucleon is excited 
from an initial single particle state a with the energy to a final one /? in the continuum with 
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the energy ep. A final expression of the inclusive {N,N'x) double differential cross section for 
the one-step process is given by 

where A is the target mass number, kc and fec(r) (c=i or /) the wave number at infinity and 
the local wave number in the initial (z) and final (/) channels, xt ixj) the distorted wave in 
the initial (final) channel, and Pt{^) {j=p or n) the nucleon density for proton {jp) or neutron 
(n) as the struck nucleon. 

The local average nucleon-nucleon (NN) scattering cross section for collision between a lead- 
ing particle, N, and a target nucleon, r, is given by 

(5(k/(r) + ka-ki(r)-k«)5(e;3-£„-'^), (2) 

where ka(k:3) is the nucleon momentum for a single-particle state with the energy ea{£p), the 
energy transfer, (da/dQ,)NT the two-nucleon scattering cross section and /^(k, r) is the Wigner 
transform of the OBDM for the hole states of proton and neutron. In Eq. (2) the sum rule of 
the WT for the hole and particle states is used: 

A(k,r)+/p(k,r) = 2, (3) 

where /p(k, r) stands for the WT for the particle states. The factor /^(ka,r)[2 — /^(k^,r)] in 

Eq. (2) represents the "probability" of the momentum state /iko, being occupied and ^k^ being 
unoccupied. The Pauli blocking effect is reflected in this factor. The difference between the 
single-particle energies in the second 5-function in the right-hand side of Eq. (2) is given by 

- ep = ^^^[k}ir) - kUr)], (4) 

using the effective mass of a nucleon, m*{r) [16]. 

The extension to higher-step processes is straightforward and the same final expressions as 
in Ref. [4] with the local average NN scattering cross sections given by Eq. (2) for successive 
collision points are also deduced for two- and three-step processes. 

3 Wigner transform of one-body density matrix using different 
theoretical models 

We introduce formulae for the Wigner transform (WT) derived in some theoretical methods, 
namely calculated with WS s.p. wave function (denoted as WS), with natural orbitals (NO) ob- 
tained from the Jastrow correlation method (JCM) which accounts for short-range correlations 
(SRC) in the case of ^'^Ca (denoted as JCM), from the approach based on the local density ap- 
proximation (LDA) including the effects of SRC (denoted as LDA-I-SRC) and from the coherent 
density fiuctuation model (denoted as CDFM). 

3.1 WT with WS 

For a single-particle potential, for instance Woods-Saxon potential, the WT for the hole states 
of proton or neutron is given by [5] 

2i -I- 1 /'°° 

41k, r) = ^ -ff - / dse-'^-^ Yl ^nlmji^ + s/2)&,(r " s/2), (5) 

nlj + ^ -^0 ^ 
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where (pnimji^) ^-P- '^^^e function for the target nucleon and the sum runs over all the 

occupied orbits nlj of protons or neutrons. 



3.2 WT with JCM 

A model independent way to define a set of single-particle wave functions and occupation prob- 
abilities uniquely from the correlated OBDM p{r, r') is to use its natural orbital representation 
(NOR) [17] 

pNOR{r, r') = naV'a(r)V'a(r')- (6) 

a 

The normalized eigenfunctions V'a(r) of p(r, r'), the so-called natural orbitals, form a complete 
orthonormal set. The associated eigenvalues called natural occupation numbers define the 
probability (0 < < 1) of the natural orbital ipa{^) occupation in the ground state Within 
the NOR the Wigner transfrm of p(r, r') can be derived as follows: 

/■oo 

/(ke,r) =5^ne/ifO(k„r) = 5^nc / dse-^^=^Vc(r + s/2)V':(r - s/2) (7) 

Wc would like to note that the summation in Eq. (7) is taken over all hole- as well as particle- 
state NO included in the calculation of the Wigner transform. In the present work we use the 
results for the natural orbitals and for the occupation numbers in ^'^Ca obtained in Ref. [18] 
within the JCM in its low-order approximation. As it has been already done in [5] for the s.p. 
wave functions, in the present work each NO Vc(r) is expanded in terms of the Gaussian-type 
basis function 

TV 

Mr) = T.''v^^M-4r% (8) 

v=l 

where N is the number of the basis functions, is given by a geometrical progression [19], = 
Kii^N The {N, Ki, kn} arc the input parameters. The expansion coefficients 
a^f^ are determined by linear fitting of the basis functions (8) to the numerical data for the 
natural orbitals. Then each partial Wigner transform is calculated in accordance with Eq. 
(2.39) in Ref. [5] and after applying Eq. (7) the total Wigner transform for the NOR of p{r, r') 
is obtained. 



3.3 WT with LDA+SRC 



It is well known that the inclusion of correlations in nuclear matter modifies the occupation 
probability predicted by the local density Fermi gas model. According to Ref. [20] one can 
introduce Wigner transform in a general way: 



f{kF{r), k) = Q{kF{r) -k) + 5f{kF{r), k), 



(9) 



where Q{kF — k) corresponds to the WT of the LEG model, while 6f{kF, k) is entirely due to 
the effects of dynamical correlations induced by the NN interaction. The local Fermi momentum 
kF{r) is related to the mass density through the relation 



kF{r) 



-7r^p(r) 



1/3 



(10) 



By definition of kF{r) one has / Sf{kF{Tc),k)d'k = 0. It was shown in [20] that for a finite 
nucleus the separation between mean-field contribution and correlation effects can be performed 
in an analogous way. For convenience more phenomenological procedure based on the results of 
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the lowest order cluster (LOG) approximation developed in [21] has been followed to evaluate 
explicitly the correlated term. Choosing a correlation function of the form 

f{r) = l-e-^''\ (11) 

the LOG gives for df{kF, k) 

df{kF{r),k) = [Y{k,8) - kdir] e{kF{r) - k) + 8 {kdirY{k,2) - [Y{k,4)]^} , (12) 

where 

c-^Y{k, fi) = -~ + / dy + sgn{k_) / e"^ dy (13) 

2k Jo Jo 



with 



The quantity 



1 //7,VV2 k ~ kF±k , , X 



...^|f/(/(.)-l)^c.r^^(|) (15) 
is the direct part of the Jastrow wound parameter. 

3.4 WT with CDFM 

The GDFM has been suggested in [22, 23, 24] as a model for studying characteristics of nuclear 
strTicturc and nuclear reactions based on the local density distribution as a variable of the theory 
and using the essential results of the infinite nuclear matter theory. The model is introduced 
using the main ansatz of the generator coordinate method for the many-body function and the 
delta-function approximation for the overlap and energy kernels of the corresponding integral 
equation for the weight function. In the GDFM the Wigner distribution function can be written 
in the form: 

/(r, k) = / |/(x)|2e(x - \r\)QikFix) - \k\)dx, (16) 
Jo 

where the weight function /(x) in the generator coordinate method is determined under the 
condition 

poo 

/ \f{x)\^dx = l. (17) 

In the case of monotonically-decrcasing density distributions {dp/dr < 0) one can obtain a 
relation of the weight function f{x) with the density distribution: 



I m2 1 Mr) 



Po{x) dr 
where 



(18) 



po{x) = SA/Attx^ (19) 

and the generator coordinate x is the radius of a sphere containing all A nucleons uniformly 
distributed in it (the so-called "flucton"). 
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4 Results of calculations and discussion 



The input data used in the calculations of Wigner transforms can be summarized as follows. 
The single-particle wave function in Eq. (5) were calculated using a Woods-Saxon potential with 
the radius parameter ro=1.27 fm and the surface diffuseness parameter a=0.67 fm, including 
the isovector term and the Coulomb potential. The set {N, ki, «;Ar}={20, 1, 0.1} is taken for the 
WS potential and for the WT with NOR. For the correlation factor in Eq. (11) we adopt the 
same value (3=1.1 fm as in Ref. [20] which is taken from the microscopic nuclear matter 
calculations [25] but reproduces also well the results for the nucleon momentum distribution 
n{k) of finite nuclei. This value leads to kdir=0.23 [Eq. (15)]. In the CDFM method as well 
as in the LDA approach with SRC we use the nucleon density distribution p(r) of the Woods- 
Saxon shape with Negele's geometrical parameters [26]. The proton and the neutron densities 
are assumed to be given by {Z/A)p(r) and (N/A)p(r), respectively. The local Fermi momentum 
kpir) was obtained from Eq. (10) using the same density distribution. 

As for the other input data for the SCDW calculations, we use basically the same ones as 
in Ref. [8]. The global optical potentials based on Dirac phenomenology of Hama et al. [27] 
for protons and Ishibashi et al. [28] for neutrons were used. For a simplicity of the numerical 
calculations we neglect the spin-orbit coupling in the distorting potentials. The nonlocality 
correction for distorted waves and WS s.p. wave functions was taken into account using the Perey 
factor with range 0.85 [29]. An effective NN interaction in terms of G matrix parameterized by 
the Melbourne group [30] was employed. The two-nucleon scattering cross section in Eq. (2) is 
calculated using this G matrix. In addition to [8], an effective mass of a target nucleon m* is 
used. We assume a simple WS form for the r-dependence of m*. The latter at r=0 was chosen 
to be m* = 0.8m with a bare mass m. 

Figure 1 shows the SCDW double differential cross sections with CDFM, LDA-I-SRC and 
WS for the ^'^C{p,p'x) reaction at 150 MeV incident and 80 MeV outgoing energies compared 
with two sets of experimental data [11, 12]. A good overall agreement with the data including 
backward angle region is obtained in the case of CDFM. The SCDW cross section calculated with 
WS drops down at backward angles, while those with LDA-f-SRC overestimates the experimental 
data. As will be discussed later, the main reason for this behavior is the different account for 
the high-momentum components of the nuclear Fermi motion. In Fig. 1 the contributions of 
individual multistep processes when CDFM is used for the same reaction are also plotted. One 
can see that proton emission via one-step process is dominant in the angular region from 20° to 
50°. Contributions of two- and three-step processes become appreciable with increasing angles. 
The role of higher-step processes has been already seen in previous calculations [4, 5]. 

Comparisons between the cross sections calculated by using the Wigner transforms obtained 
with all theoretical methods considered and the measured cross sections for {p,p'x) on ^'^Ca at 
186 MeV and 392 MeV incident energies are given in Fig. 2 and Fig. 3, respectively. Both 
figures confirm the fact that the cross sections at large angles depend strongly on the theoretical 
model used to describe the nuclear states. Fig. 2 shows a fair agreement of the CDFM and JCM 
results with the data for emission energy of 110 MeV and to less extent for energy of 130 MeV. 
At the same time, for emission energy of 150 MeV the WS result is able to reproduce the angular 
distribution at backward angles. Similar behavior can be seen from Fig. 3, where the double 
differential cross sections calculated with CDFM and that with JCM are in a good agreement 
with the experimental data. However, measTirements over an extended range of backward angles 
can test the various methods and, hence, the role of NN correlations. 

The improved description of the backward proton emission with the SCDW model by the 
use of realistic models for nuclear states can be explained by the difference in the momentum 
distribution of the target nucleons. Figure 4 shows the momentum distributions n{k) of ^^C and 
^°Ca with CDFM, LDA-^SRC, WS and JCM. One sees that n{k) calculated with aU methods 
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Figure 1: Comparison of SCDW calculations with measured double differential cross section for 
the reaction ^^C(p, p'x) at 150 MeV incident and 80 MeV emission energies (solid circles from 
[11] and solid triangles from [12]). The calculations with CDFM (sohd line), LDA+SRC (dotted 
line) and WS (dash-dotted) line are shown in the left panel. The CDFM calculations of the cross 
sections of one- (dotted line) , two- (dash-dotted line) and three-step (dashed line) processes are 
shown in the right panel, where the solid line corresponds to their sum. 

but not with the WS contain high- momentum components at momenta larger than 1.5 fm~^. 
The LDA-I-SRC [20] and Jastrow method [18] account for short-range correlations, while the 
CDFM [22]takes into account mostly long-range correlations (LRC). The larger high-momentum 
component due to the stronger SRC when using LDA-I-SRC causes the larger cross sections at 
backward angles. Although n{k) calculated in this method fits well the experimental y-scaling 
data for ^^C [31] at large momenta, the corresponding cross sections overestimate the data. 
Here we would like to note that the n{k) of ^°Ca obtained after integration of the total WT 
/(k, r) with JCM over r is not given in Fig. 4 and it underestimates the one obtained in [18] at 
k > 2 fm~^. This is particularly due to the truncation of the particle-state NO's included in the 
calculation of the WT of ^Ca. On the other side, due to the lack of enough high-momentum 
components the angular distributions calculated with WS underestimate the data at backward 
angles. The nucleon momentum distributions obtained with CDFM for both nuclei presented 
in Fig. 4 also exhibit high-momentum components, though smaller than the ones corresponding 
to methods accounting for SRC. Fig. 4 clearly illustrates the strong dependence of backward 
proton production from {p,p'x) reactions on the momentum distribution of target nucleons. 

In Fig. 5 the SCDW cross sections for ^^Zv{p,p'x) at 160 MeV and emission energy of 120 
MeV calculated with the CDFM, LDA-I-SRC and WS models are compared. Again a very good 
agreement with the experimental data over the forward, intermediate and backward angles is 
achieved when using the Wigner transform obtained from CDFM. In Fig. 6 the contributions of 
individual multistep processes for the same reaction and corresponding to the same models are 
plotted to show their variation with the scattering angle. It is seen that for this relatively large 
emission energy the one-step cross section calculated with LDA-j-SRC exceeds the experimental 
cross section for large angles. On the other hand, the one-step cross section obtained with 
WS drops at angles less than 140° underestimating the experimental cross section at backward 
angles. 
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Figure 2: Comparison of SCDW calculations with measured double differential cross section for 
the reaction ^^C&{p,p'x) at 186 MeV incident and 110 MeV (a), 130 MeV (b) and 150 MeV (c) 
emission energies [13]. In each panel the calculations with CDFM (solid line), JCM (dash-double 
dotted line), LDA+SRC (dotted line) and WS (dash-dotted line) are shown. 

From the analyses of high energy heavy ion collision and deep inelastic lepton-nucleus scat- 
tering, the following form of the nucleon momentum distribution has been proposed [10]: 

n{k) = N[exp{-kypl) + £0 exp(-A;Vgo) + exp{-k'^ /qj)], (20) 

where iV is a normalization factor and the values of the parameters are determined to be po = 
^J\kF, eq ^ 0.03, ^ v/3po, £1 = 0.003 and gi ^ 0.5 GeV/c. The first term of Eq. (20) 
corresponds to single-particle contribution to n{k) in the mean-field approximation (MFA). It 
was noticed in [10] that the second term in Eq. (20) is due to the long-range correlations (L- 
type high momentum component) while the third term was supposed to be due to short-range 
correlations (S-type high momentum component). Moreover, as it was also noticed in Ref. 
[10] only the high-momentum component generated by the Hartree-Fock correlation should be 
included in the single scattering calculation. In order to reveal better the role of different NN 
correlations for the correct description of the backscattered proton spectra, we make in Fig. 7 
a comparison between the nucleon momentum distributions calculated with CDFM and JCM 
and those calculated by using Eq. (20) in the case of ^^C and ^'^Ca nuclei. Particularly, our 
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Figure 3: Comparison of SCDW calculations with measured double differential cross section for 
the reaction '^^Ca.{p,p'x) at 392 MeV incident and 100 MeV (a), 200 MeV (b) and 300 MeV 
(c) emission energies [14]. In each panel the calculations with CDFM (dotted line), JCM (solid 
line), LDA+SRC (dash-dotted line) and WS (dashed line) are shown. 

interest is in the comparison of the high-momentum components which the different methods 
contain. It turns out from Fig. 7 that for k > 1.5 fm^^ the curves corresponding to the inclusion 
only of the L-type term and CDFM have similar behavior (especially for ^^C), while the results 
when including only the S-type term and Jastrow method are close to each other. It is seen also 
from Fig. 7 that the contribution of the S-type term dominates in the total n{k) calculated by 
Eq. (20). Therefore, we suggest that the long-range correlations play an important role to the 
problem of backward proton production. This is confirmed by the good agreement obtained in 
[9] and by our results obtained with the use of Wigner transform with CDFM which lead to 
best description of the experimental cross sections over wide target muss number and emission 
energy regions. 

5 Summary 

The results of the present work can be summarized as follows: 

i) Double-differential cross sections of ^'^C{p,p'x), '^^Ca{p,p'x) and ^^Zr(j),p'x) reactions are 
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Figure 4: Comparison of nucleon momentum distributions n{k) (normalized to unity) of ^^C 
and ^°Ca calculated with CDFM (dotted line), LDA+SRC (solid line), WS (dashed hue) and 
JCM [18] (dash-dotted line). The experimental data for n{k) of ^^C are taken from [31]. 

calculated, using SCDW model and Wigner transform of the OBDM obtained with different 
approaches, in the incident energy range of 150-392 MeV. A good overall agreement with the 
experimental data is achieved including the backward angles region. These analyses show that 
the present SCDW model can describe successfully the MSD processes of {p,p'x) reactions over 
the wide range of target mass number and incident energies. 

ii) It is found that the high-momentum components of the nucleon momentum distribution 
are most responsible for describing the angular distributions of the backscattered protons. Since 
CDFM leads to the best agreement with the experimental cross sections at wide emission energy 
region, it is demonstrated that the long-range correlations which are related to the collective 
nucleon motion are rather important than the short-range ones generated in the other correlation 
methods. 

A more detailed study on the backward proton emission in the framework of the CDFM is 
in progress. 
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